Numerical simulation of fracture mode transition in ductile plates

Fracture mode of ductile solids can vary depending on the history of stress state the material experienced. For example, ductile plates under remote in-plane loading are often found to rupture in mode I or mixed mode I/III. The distinct crack patterns are observed in many different metals and alloys, but until now the underlying physical principles, though highly debated, remain unresolved. Here we show that the existing theories are not capable of capturing the mixed mode I/III due to a missing ingredient in the constitutive equations. We introduce an azimuthal dependent fracture envelope and illustrate that two competing fracture mechanisms, governed by the pressure and the Lode angle of the stress tensor, respectively, exist ahead of the crack tip. Using the continuum damage plasticity model, we demonstrate that the distinctive features of the two crack propagation modes in ductile plates can be reproduced using three dimensional finite element simulations. The magnitude of the tunneling effect and the apparent crack growth resistance are calculated and agree with experimental observations. The finite element mesh size dependences of the fracture mode and the apparent crack growth resistance are also